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Abstract 

We have written the Exclusive Hadronic Monte Carlo Event (Ex- 
HuME) generator. ExHuME is based around the perturbative QCD cal- 
culation of Khoze, Martin and Ryskin of the process pp — > p + X + p, 
where X is a centrally produced colour singlet system. 

Introduction 

Central exclusive production events are those in which the colliding protons re- 
main intact after an interaction and all of the energy transferred from the pro- 
tons goes into a central colour singlet system. Unlike inclusive double pomeron 
exchange, in which pomeron remnants are present, there is no radiation other 
than that from the central system. Furthermore, in the limit that the transverse 
momentum of each of the outgoing protons vanishes, no angular momentum is 
exchanged between the proton lines and the quantum numbers of the central 
system are constrained to be Jz = with even parity. The leading order dia- 
gram for the central exclusive production of system X is shown in figure ^ 
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Figure 1: The Exclusive Diffractive Process 

The invariant mass of the central system can be determined if proton de- 
tectors mounted very forward of the central detector are used to measure the 
longitudinal momentum losses, xi,2, and momentum transfers squared, ti t 2- 

ExHuME factorises central exclusive production as indicated by the dashed 
line in figure^ We use gluon fusion cross sections - taking the colour singlet and 
spin zero projections of the amplitude pQ - but do not fix the collision energy. 
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The luminosity, C, for the fusing gluons is dependent on both the invariant mass 
(V§) and rapidity, (y = ^ log(a;i/x2)) of the central system QQ[21EI|: 



=eb(tl+t2) (i / ^ /s(Q± ' M ' xi > ^^x, i", * 2 , 4))' 



dsdydtidt 2 



(1) 

The forward limit implies that feu_ ~ — fc2± — — Q±, which simplifies the integral 
to the form in equation QJ. The f g are the unintegrated skewed gluon distri- 
butions (which can be calculated from the integrated gluons plj). We choose 
fx = 0.618V 7 ! for the hard scale and 6 = 4 GeV -2 pp. The x' are the longi- 
tudinal momentum fractions of the gluon that is not connected to the central 
system, x' x = x' 2 = 0{Q 2 L / s) for symmetric proton beams of energy yfs/2. The 
integral in equation Q converges in the infra-red due to the presence of a Su- 
dakov factor that imposes the requirement that there be no radiation emitted 
between the scales of Q± and fx. The leading term in ki± ■ k 2 ± in the integrand 
of equation is Q\ and has no dependence on the azimuthal angle cj> between 
the outgoing proton directions. For complete details of the luminosity calcula- 
tion we refer to the calculation of V. A. Khoze et al [2j EJ Hj ; upon which 
ExHuME is based. 

A survival factor, S 2 , is also required to ensure that there are no additional 
interactions between the proton lines. We use a constant value for S 2 that 
is fixed for a given collision energy, although we anticipate that further work 
will lead to a parameterisation of the survival factor in terms of t\, t 2 and <fi 
[1]. We note that the 4> dependence in that model is approximately flat for 
values of P3j_, p4± < 200MeV. The default differential luminosity curves used 
by ExHuME are shown in figure EJ 

In this initial release of ExHuME we provide the following sub-processes: 
gg — > H, gg — > QQ and gg — > gg, where H is a Standard Model Higgs and 
Q is a massive quark. We require a width for resonance production and use 
the Higgs line shape described in j^] with the Higgs width and branching ratios 
calculated using the program Hdecay [Jj 
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Figure 2: The default luminosity used in ExHuME at collider energies of 1.96 
TeV (dashed line) and 14 TeV (solid line), a) shows the luminosity vs. the 
invariant mass of the fusing gluons at zero rapidity, b) shows the rapidity 
dependence of the luminosity when the invariant mass of the fusing gluons is 
fixed at 120 GeV. A fixed survival factor of 0.03 (14 TeV) or 0.045 (1.96 TeV) 
has been included and the luminosity has been integrated over the p± of the 
outgoing protons. 



Design 

ExHuME is written in a modular way using C++. There are two main classes 
that are needed to generate central exclusive events. The first is a CrossSection 
class that calculates the differential luminosity, the gluon fusion sub-process and 
the kinematics of any outgoing particles. The second class needed is an Event 
class that generates the events. In addition there is a Weight class that provides 
random numbers distributed according to any given function and from which 
the Event class inherits. 
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The CrossSection Class 



The CrossSection class exploits the factorisation of figure^] CrossSection is 
an abstract base class containing the calculation of the effective luminosity of 
the gluon-gluon collision with a virtual method for the gluon fusion sub-process. 
Complete processes are created by inheriting from the CrossSection and ex- 
plicitly defining a sub-process. This makes it relatively simple to implement 
new sub-processes in addition to the currently implemented Standard Model 
Higgs, gg and QQ sub-processes. Instructions for creating a new sub-process 
are posted at jH] and users are invited to submit any new process to the authors 
for inclusion in future updates to ExHuME. 

CrossSection Methods 

The constructor for the CrossSection class is 

CrossSection(int , char**), 

which allows the user to pass a card file from the command line that changes the 
default values used in the luminosity calculation (see Appendix EJ. The con- 
structors for the derived classes are Higgs (int, char**), QQ(int, char**), 
GG(int , char**) orDummy(int, char**). 

The invariant mass and rapidity of the central system and the momentum 
transfers and azimuthal angles, </>i,2, of the outgoing protons can be set by the 
CrossSection method 

void SetKinematics ( 

const double femass, const double ferapidity, 
const double &tl, const double &t2, 
const double ftphil, const double &phi2);. 

The method 

double Dif f erential () 

then returns the differential cross section. The 4- vectors of the outgoing protons 
can be accessed by the methods 

HepLorentzVector GetProtonlO and 
HepLorentzVector GetProton2 () . 

The sub-process level information can be extracted by the method 
std: : vector<Particle> GetPartonsO 

which returns the outgoing particles from the gluon fusion matrix element 
prior to decay, parton showering and hadronisation. The Particle class that 
GetPartonsO returns contains the particle momentum, p, the PDG ID code, 
id and vertex, vtx. It is also possible to access the kinematics of the cross 
section by the following methods: 
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double GetRoot_s(); 

Returns the invariant mass of the colliding beams, 
double GetsHatO; 

Returns the invariant mass squared of the central system, 
double GetSqrtsHat () ; 

Returns the invariant mass of the central system. 

double GetxlO; 
Gets the value of x±. 

double Getx2(); 
Gets the value of X2- 

double GettlO; 
Gets the value of t\ . 

double Gett2(); 
Gets the value of t% ■ 

double GetPhilO; 

Gets the value of <p\ , the azimuthal angle between the outgoing proton 1 and 
the x direction. 

double GetPhi2(); 

Gets the value of </>2, the azimuthal angle between the outgoing proton 2 and 
the x direction. 

double GetEtaO; 

Gets the rapidity, y, of the central system. 
Hadronisation is performed via the method 
void HadroniseO 

which places the event record into an external hepevt common block that can 
be written to file. 

There are a number of sub-process specific methods. The Higgs decay type 
can be set by the method 

SetHiggsDecay (const int&) 

with the PDG code of the decay products as the argument. Similarly the quark 
type in qq can be set by the method 

SetQuarkType (const int&) 

which defaults to bb production. In addition, the GG and QQ classes contain 
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the method 

SetThetaMin( const doublefe) 

which sets the minimum (and maximum) polar angle of an outgoing parton 
relative to the beamline in the rest frame of the central system. The default 
value is cos 9 = 0.95. 

The Event Class 

The Event class generates events, calculates the total cross section and reports 
the efficiency with which events were generated once event generation has fin- 
ished. Event contains a pointer to a CrossSection and the user specifies which 
sub-process is to be generated in the constructor for the Event. The Event 
class also inherits from the Weight class. The numerical weighting algorithm 
is initialised to return the mass distribution of the differential cross section at 
a rapidity zero and is effective even for a narrow resonance such as the Higgs. 
The variables t\ and t 2 are distributed according to e b ^ tl+t2 ^ whilst </>i, 4>2 and 
y are uniformly distributed. 

Event Methods 

The event is defined by calling the constructor 
Event (CrossSectionfe P, const unsigned int R) , 

where P is a CrossSection with the sub-process defined and R is a random 
number seed. There are a number of methods that can be used to set the kine- 
matic ranges of the parameters used to define the event: 

void SetxlMax(const doublefe); 
Sets the upper limit of Xi. 

void Setx2Max( const doubleft); 
Sets the upper limit of x-i . 

void Sett lMax( const doubleft); 
Sets the upper limit of t\ (t < 0). 

void SettlMin(const doublefe); 
Sets the lower limit of ii . 

void Sett2Max(const doubleft); 
Sets the upper limit of t 2 - 

void Sett2Min(const doublefe); 
Sets the lower limit of t 2 ■ 

void SetMassRange (const double & minimum, const double & maximum); 
Sets the lower and upper limits of the mass range respectively. 
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The last method that is used before event generation is 
void SetParameterSpace () , 

which must be called in order to initialise the event generation. Individual 
events are generated by the method 

void Generate (). 

Finally, the total cross section from the events generated can be calculated 
by the method 

double CrossSectionCalculationO 
and the efficiency via 
double GetEf f iciency () . 

Using ExHuME 
Installing ExHuME 

The ExHuME source code is available from jH] or on request from the authors. 
In its standard form ExHuME must be linked at compilation to Pythia Pj, 
CLHEP [12] and either LHAPDF p] or the CERN PDFLIB [TJ]. It would 
also be possible to modify ExHuME to use Herwig JOj instead of Pythia for 
the hadronisation or to use a stand alone PDF instead of either LHAPDF or 
PDFLIB. By default ExHuME sets the location of the directory containing the 
grid or parameter files for LHAPDF to be wherever the program is executed 
from. A symbolic link should be created to wherever the grid and parameter files 
actually reside. For further information please see the respective documentation 
for each of these programs . 

Example Main Program 

In this section we demonstrate a simple main program that generates 5000 
H — > WW* events for a Higgs with the default mass of 120 GeV. We also show 
how to extract simple information from the hepevt record. 

#include <iostream> 

The following headers are for ExHuME: 

#include "Event. h" 
#include "Higgs. h" 

int main(int argc, char** argv){ 



7 



Declare a new Higgs CrossSection: 

Exhume: : Higgs higgs (argc , argv) ; 
and set the Higgs to decay only to W bosons: 

higgs. SetHiggsDecay (24) ; // 24 is the PDG code for W 

Declare an event with the Higgs as the cross section and a random number 
seed of 1111. 

Exhume: : Event HiggsEvent (higgs , 1111) ; 

The allowed range of gluon fusion invariant masses must be set. As long as the 
range is much bigger than the width of the resonance the results will not be sen- 
sitive to the range chosen. This is, of course, not the case whenever the central 
system does not have a narrow width, for example in gg and QQ production. 
The efficiency of event generation will drop if a large mass range is chosen for 
a narrow resonance such as the Higgs, so for this example where the width is 
0.0036 GeV we set the mass range to be between 115 GeV and 125 GeV. 

HiggsEvent . SetMassRange(115, 125) ; 

This must be called before event generation can begin: 

HiggsEvent . SetParameterSpace () ; 

double xl; 
int Nobj ; 

std: : vector<Exhume : :Particle> Higgslnfo; 
for ( int i = ; i < 5000 ; i++ M 
The next line generates a single event: 

HiggsEvent . Generate () ; 
Access the longitudinal momentum loss of proton 1 for this event: 

xl = process . Getxl () ; 
and the information about the Higgs: 

Higgslnfo = process . GetPartons () ; 
Get the number of particles in the hepevt common block: 
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Nobj = hepevt_ .nhep; 

> 

std: : cout<<" Cross section = " 

<<HiggsEvent . CrossSectionCalculationO <<std: : endl ; 
std::cout<<" Efficiency of event generation = " 
<<HiggsEvent . GetEf f iciency () <<std: : endl ; 
return (0) ; 

} 

The program allows (but does not demand) a card file to be given on the com- 
mand line that overrides the default parameters. Such a card file could look like 

HiggsMass 140 
TopMass 180 

which would be appropriate for investigating the effects of varying the higgs 
and top masses. 
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A The ExHuME Control Parameters 

ExHuME can be controlled by passing a card file from the command line that 
can contain the control parameters given in table The collider can be set 
to the LHC (1), the Tevatron (0) or neither (—1). Choosing a collider sets the 
proton collision energy, yfs, the survival factor S 2 and R g , which accounts for 
the skewed effect in the un-integrated gluons. Choosing neither means that the 
user must set these parameters. The PDF values are the PDF set numbers 
accepted by the LHAPDF library. Freeze is a scale below which a s is frozen 
and Q± m in is the lower bound on the integral in equation (QJ - the integral is 
formally convergent, but the necessary freezing of a s requires a lower bound to 
be introduced. 
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Parameter 


Name 


Type 


Default 


a 


AlphaEW 


double 


0.0072974 


M w 


WMass 


double 


80.33 


M z 


ZMass 


double 


91.127 


M H 


HiggsMass 


double 


120.0 


M t 


TopMass 


double 


175.0 


M b 


BottomMass 


double 


4.6 


M c 


CharmMass 


double 


1.42 


M s 


StrangeMass 


double 


0.19 


M T 


TauMass 


double 


1.77 




MuonMass 


double 


0.1057 


V 


HiggsVev 


double 


246.0 


Q 2 , . 

™ 1mm 


MinQt2 


double 


0.64 


Aqcd (MeV) 


LambdaQCD 


double 


80 (MeV) 




Freeze 


double 


Q-Lmin 


b 


B 


double 


4.0 


collider 


FNAL_or_LHC 


int 


1 


s 


s 


double 


1.96 x 10 8 


i 

S2 


root_s 


double 


14000.0 




Rg 


double 


1.2 


5 2 


Survive 


double 


0.03 


PDF 


PDF 


int 


20250 



Table 1: ExHuME control parameters 
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